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The Majorana code is an example of a stabilizer code where the quantum information is stored in 
a system supporting well-separated Majorana Bound States (MBSs). We focus on one-dimensional 
realizations of the Majorana code, as well as networks of such structures, and investigate their 
lifetime when coupled to a parity-preserving thermal environment. We apply the Davies prescription, 
a standard method that describes the basic aspects of a thermal environment, and derive a master 
equation in the Born-Markov limit. We first focus on a single wire with immobile MBSs and 
perform error correction to annihilate thermal excitations. In the high-temperature limit, we show 
both analytically and numerically that the lifetime of the Majorana qubit grows logarithmically 
with the size of the wire. We then study a trijunction with four MBSs when braiding is executed. 
We study the occurrence of dangerous error processes that prevent the lifetime of the Majorana 
code from growing with the size of the trijunction. The origin of the dangerous processes is the 
braiding itself, which separates pairs of excitations and renders the noise nonlocal; these processes 
arise from the basic constraints of moving MBSs in ID structures. We confirm our predictions with 
Monte Carlo simulations in the low-temperature regime, i.e. the regime of practical relevance. Our 
results put a restriction on the degree of self-correction of this particular ID topological quantum 
computing architecture. 
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I. INTRODUCTION 


Topological quantum computation (TQC) describes 
the general idea of storing and processing quantum infor¬ 
mation in topological states of matter.^ The most ap¬ 
pealing aspects of TQC reside in the intrinsic protection 
of the ground-state subspace against local (static) per¬ 
turbations; topological ground states are thus viewed as 
a good place to hide quantum information. Furthermore, 
quantum gates are executed by performing highly non¬ 
local operations that consist in the exchange (or braid¬ 
ing) of quasi-particles in the form of non-abelian anyons. 
While it is difficult for the environment to induce such 
exchanges, an external observer is able to do it by adia- 
batically controlling the parameters of the system. Also, 
the applied quantum gates depend only on the topology 
of the exchange and are thus insensitive to local imper¬ 
fections. 


In the last decade, it has appeared that Ising anyons 
are the non-abelian particles most likely to occur in phys¬ 
ical systems in the laboratory, see Ref. [2] and references 
therein. Although their braiding statistics is not rich 
enough to generate a universal set of gates, they allow 
the implementation of the Clifford group in a topologi¬ 
cally protected fashion and are thus of strong interest for 
quantum computation. In this context, the so-called Ki- 
taev wir^^has recently attracted tremendous attention. 
In fact, unpaired Majorana modes appear in this model 
and, when braided in a network of one-dimensional wire^ 
they behave as IsIm anyons .1^ Considerable theoreticaP^^ 
and experimentaP^Hi^ efforts have been invested to inves¬ 
tigate semiconducting hybrid structures that could real¬ 


ize the Kitaev wire. 

Although Majorana qubits exhibit many favorable 
properties, more and more studies have focused on the 
fragility of such topological qubits. In particular, several 
sources of noise that limit the applicability of such setups 
have been reported.^^^^^ 

In this paper we start from a microscopic model and 
study the lifetime of the Majorana code, see Rnfs. fTbllSTl 
and [3^ as well as Sec. |IIB| for a definition, when coupled 
to a parity-preserving thermal environment. We apply 
the Davies prescription to derive a Born-Markov master 
equation. We first focus on a single wire with immo¬ 
bile Majorana Bound States (MBSs) and discuss how to 
perform error correction to counteract the effect of the 
environment. We demonstrate in the high-temperature 
limit, both analytically and with Monte Carlo methods, 
that the lifetime grows logarithmically with the system 
size. This result is not unexpected as similar behavior 
was observed by Bravyi and Ko enig for a closed system 
with Hamiltonian perturbations.^ As a next step, we 
study a trijunction with moving MBSs. Our main result 
is the investigation of details of dangerous error processes 
that prevent the lifetime of the system from increasing 
with the system size. The origin of dangerous errors is 
the braiding itself that renders a local error source highly 
non local by dragging excitations across the trijunction. 
In particular, we demonstrate that performing error cor¬ 
rection at the end of the braid does not allow the danger¬ 
ous errors to be cured. We confirm our predictions with 
a Monte Carlo simulation. Our work is an extension of 
Ref. [32| Here, we present additional physical results as 
well as the technical details leading to the main results 
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of Ref. |32j 

In the context of a full quantum computing protocol, 
where several braids are executed, our results imply that 
error correction at the end of all the braids, i.e. purely 
passive, is not enough. Our results show also that a more 
active scheme, in which error correction is executed at the 
end of each braid, is also too weak to solve the problem 
of dangerous errors. We are led to the view that error 
correction will only be successful if it is fully active, i.e., 
where several error correction steps are executed dur¬ 
ing each braid, to counteract the decoherence effects of 
dangerous errors. Therefore, our work brings additional 
evidence that even in non-abelian topological codes ac¬ 
tive error correction, in the same sense as for ordinary 
quantum error correction codes, is necessary.l^^HllI 

The paper is organized as follows. In Sec.|n|we present 
the main aspects of a single Kitaev wire that carries 
MBSs at the junction between topological and nontopo- 
logical segments. In Sec. II A| we introduce a box repre¬ 
sentation of the wire that turns out to be useful to under¬ 
stand the phenomenology of the wire as well as the way 
we simulate it. In Sec. |IIB| we define the Majorana code, 
i.e. a stabilizer code that encodes a logical qubit in the 
ground-state subspace of the Kitaev wire. In Sec. |II C 


we 


define string operators that create, annihilate, and move 
excitations in the wire. The string operators give us a 
rigorous way to perform error correction. In Sec. IIP we 


study the coupling between the Kitaev wire and a bosonic 
bath. We follow the Davies prescription and derive a 
Markovian master equation in Sec. nil In Sec. m we 
focus on the lifetime of the single wire Majorana code at 
high temperatures. We derive an analytical formula for 
the lifetime in Sec. HIE H and confirm it with Monte Carlo 
simulations in Sec. lIIEH In Sec. Hill we introduce the tri¬ 


junction setup used to braid MBSs and in Sec. Ill A we 
show how the logical qubit is encoded in four well sepa¬ 
rated MBSs. In Sec. [IV| we study in detail the unitary 
evolution arising when MBSs are moved. In particular, 
we focus on the behavior of excitations. In Sec. |IV A[ we 
present a rigorous definition of what adiabaticity means 
in our study. In Sec.|V|we show how the master equation 
for the time-independent Kitaev wire generalizes to the 
time-dependent trijunction setup in the adiabatic limit. 
In Sec. |VI| we present the algorithm we use to perform 
error correction in the trijunction. Section [Vl A| contains 
our main results; we identify dangerous error processes 
that prevent the lifetime of the trijunction logical qubit 
from increasing with system size when braiding is exe¬ 
cuted. Einally we confirm our predictions with Monte 
Carlo simulations in Sec. 1 ^ The Appendices contain 
additional information and details about the derivations. 


II. SINGLE WIRE 

We review here the basic aspects of the physical model 
considered here and already exposed in our previous 
work, see Ref. [32l 


We start our study with a single wire of size L sup¬ 
porting immobile MBSs. The wire Hamiltonian i^^ 

L L-l 

i=i 

L-l 

+ -f- , (1) 

where aj and aj are fermionic creation and annihila¬ 
tion operators at site j. The first term describes a 
site-dependent chemical potential jij ^ 0. The second 
and third terms describe respectively nearest-neighbor 
hopping with t > 0 and superconducting pairing with 
A = |A|e*^. 

To understand the physics of in simple terms, it 
is useful to go to a representation in terms of Majorana 
operators, aj = —( 72^-1 +i 72 j) with { 7 i, 7 j} = 25ij 
and 7 J = 7 ^. Eor the case t = |A| and jj^j = OVj, we 
obtain the simplified expression 


L-l 

= -|A| i72i+i72i • (2) 

i=i 

The first Majorana mode 71 as well as the last Majo¬ 
rana mode 72 l are decoupled from the rest of the chain 
and [i 7 ^^, 7 i] = [^w^, 72 l] = 0- This allows one to 
define a zero-energy delocalized fermionic mode with an¬ 
nihilation operator 


= ^(71 +n 2 L) • (3) 

Using the eigenmode operators dj = 5 ( 72 ^ +^ 2 ^+ 1 ), the 
wire Hamiltonian takes the fully diagonal form 


L-l 

= (4) 

j=0 

where eo = 0 and ej = |A| for j = — 1 . As 

originally proposed by Kitaev,!^ it is tempting to encode 
a qubit in the ground-state subspace of i^w The reason 
is that local (static) perturbations lead to a ground-state 
splitting exponentially small in L. Therefore, the de¬ 
coherence induced by such undesirable splitting can be 
exponentially suppressed by increasing a parameter that 
is easy to control, namely the size of the wire. 

Away from the limit jij = 0, MBSs localized at the 
end of the chain persist as long as |/ij| ^ 2 t; we call 
this the topological phase. However, when iij ^ 0 the 
MBSs are not localized anymore on a single site but have 
support in the bulk of the chain; the amplitude of the 
MBS wave function decreases exponentially away from 
the end sites. Eor \fi\ ^ 2t the localized modes disappear; 
this characterizes the nontopological phase. Deep in the 
nontopological phase, with \ jaj \ ^ |A| = 0, the Majorana 
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Hamiltonian approaches 

= ( 5 ) 

i=i 

In Eq. (§, the Majorana modes are paired in a shifted 
way as compared to the topological case, see Eq. We 
present a pictorial representation of these two different 
pairings in Eig. 



FIG. 1. a) Pictorial representation of the Kitaev wire with 
topological (black) and nontopological (gray) segments. The 
large black dots describe the fermionic sites, while the line in- 
between describe hopping and superconducting pairing. The 
smaller dots below represent the Majorana modes whose pair¬ 
ing are depicted by lines connecting the dots. Pairings in the 
topological segment are shifted as compared to pairings in 
the nontopological segment. A possible pattern of '0- and 
'^'-excitations is shown, b) Value of the chemical potentials 
fjij corresponding to the situation in a). The chemical poten¬ 
tials in the nontopological segment have a gradient in order 
to localize the -^'-excitations. 

Having in mind the Majorana pairing pattern in the 
topological and non topological segments, it is straight¬ 
forward to see that MBSs appear at the junctions be¬ 
tween topological and nontopological segments of the 
wire, see Eig. By varying the chemical potential, 
one can increase or decrease the size of the nontopologi¬ 
cal segments and thus move the position of the localized 
MBSsP This idea will be used later in Sec. Illll to braid 
the MBSs. 


A. Box representation of the wire 

It is useful to use a box representation of the wire to 
understand the phenomenology of the model and the way 
we will simulate it. In Eig. we present the details of 
our representation; if the wire contains L sites, its box 
representation contains L + 1 boxes. It is worth pointing 
out that a box is used to represent either a Majorana 
mode or a fermionic mode. While the size of the boxes 
vary in Eig. for clarity, the size of each box has no 





nontopological 
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topological 


FIG. 2. Box representation of the Kitaev wire. A wire of 
length L is represented by L+1 boxes. A box represents either 
a full fermionic mode or a Majorana mode. Here and below, 
the boxes in the non topological sections are colored gray. 
The blue links carry an excitation and we draw a -0 inside 
the corresponding box. A -0 inside an MBS describes a flip 
of the logical parity So. a) The wire is completely topological 
with two MBSs at the ends and four excitations, b) The wire 
carries a topological and a nontopological segment, with two 
MBSs at the boundary of the topological segment. Here the 
topological segment hosts two excitations. 


meaning. In the rest of this work, all the boxes will have 
the same size. 


B. Majorana Code 

Following the approach of our previous work Ref. | 32 l it 
is convenient to take an information-theoretical approach 
to the encoding of log ical qubits into the ground states of 
the Kitaev wire.^^^^ The g round-state subspace of 
forms a stabilizer code^^^^with stabilizer operators Sj = 
^72j+i72j he., the terms in the Hamiltonian Eq. As 
usual for a stabilizer quantum error correcting code, two 
logical qubit states |0) and |1) have the property 

5 ,-| 0 ) = |0) and 5 ,-|I) = |I). (6) 

The Majorana code can then be interpret ed as a one¬ 
dimensional version of Kitaev’s toric code.^^^^ Excita¬ 
tions above the ground states are localized and defined 
through the conditions Sj = — 1 . These excitations are 
denoted as quasi-particles - 0 . We represent the code with 
boxes where the first and last boxes host the MBSs, while 
the other sites support either vacuum {Sj = +1) or a -0 
{Sj = - 1 ), see Fig.[^. We represent a flip of the logical 
parity Sq = ijij2L —So by drawing a -0 inside the left 
or right MBS. A -0 inside an MBS does not correspond 
to a real excitation since it does not cost any energy to 
be created, rather it signifies that the logical qubit has 
been flipped. A -0 inside the left MBS is proportional to 
a X Pauli, a ^ inside the right MBS is proportional to a 
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Y Pauli, and a inside both MBSs is proportional to a 
Z Pauli. 


C. String Operators, Fusion, and Error Correction 

We solely consider parity-conserving perturbations and 
particles are thus always created in pairs. Pairs of 
excitations are generated by string operators; a string 
operator creating excitations Sj = —1 and S'/. = — 1 reads 

Sjk = 72i+l72i+2 • • • 72/c , (7) 

see Fig. We define the weight of a string operator Sj k 
as \j — It is then clear that 


Z oc Sol = 7i72 * • •72L (8) 

maps the ground-state subspace into itself creating a pj 
inside the left MBS and a pj inside the right MBS. Since 
X and Y logicals are generated by an odd number of 
pj particles, they cannot be implemented in a parity¬ 
preserving scenario with immobile MBSs. As we will see, 
when MBSs are braided the situation changes drastically. 


O 0-0 O ~0 0~0 O ~0 0-0 0-0 0 -^ 00 00 0-0 00 0 

V__ J 


So 3 = 717273747576 


i 


Soo — 7l37l47l57l67l77l8 


to conserve the information stored in the ground states of 
one needs to define a protocol for error correction 
based on the knowledge of the positions of in the bulk 
(recalling that inside an MBS is invisible), the so-called 
error syndrome. If a pair of t/^’s is created in the bulk 
of the chain and not annihilated, then one particle can 
diffuse to the left end, while the second one diffuses to 
the right end. The operation performed on the logical 
qubit is then proportional to Z. 

The goal of error correction is to counteract the effect 
of the environment by finding a procedure that annihi¬ 
lates the excitations in a definite manner such that the 
stored quantum information is retrieved. Since it is rea¬ 
sonable to assume that nearby t/^’s originate from the 
same error event (as is the case at small times), we anni¬ 
hilate them following a Minimal Weight Perfect Match¬ 
ing (MWPM) algorithm for the single wire. In one di¬ 
mension there are only two possibilities to perform such 
pairingsOne of them will lead to a successful logical 
qubit recovery, while the second one will introduce a Z 
error, see Fig. Which of the two schemes is chosen 
depends on the total number of moves to be applied. We 
choose the scheme with minimal weight. Note that we 
will eventually use a different algorithm when we study 


the trijunction, see Sec. VI 
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FIG. 3. Excitations created by application of string operators 
Sj k defined in Eq. 0- 

String operators give us a way to fuse excitations 
and thus to perform error correction. Two ip particles 
Sa = —I and Sb = —1 are fused by applying Sab- The 
effect is to bring back the system into its ground state 
by annihilating the quasi-particles. Similarly, a ip par¬ 
ticle Sa = —I can be fused to the left (right) MBS by 
applying 5oo {SaL)- 

In light of the above considerations, it is clear that the 
phenomenology of the M ajor ana wire is the same as the 
Ising anyon model,I^IE21 

'0X'0 = 1, a X pj = pj ^ and a x a = 1 ip , (9) 

where a is the standard label for an Ising anyon and 1 for 
vacuum. Here a particles are identified with the MBSs. 
The second Eq. 0 indicates that, as we have seen, a ip 
inside an MBS is invisible to an external observer. Also it 
is clear that two MBSs give us a two-dimensional Hilbert 
space, cr X cr = 1 as we have seen 1 corresponds to an 
empty delocalized mode with d\do =0 and to a filled 
delocalized mode d\do = 1. 

In the following, we assume that the wire is in contact 
with a thermal bath that generates excitations. In order 
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weight 9 

EIG. 4. Two examples of the Minimal Weight Perfect Match¬ 
ing (MWPM) algorithm. The total weight of the algorithm is 
defined as the sum of the weights of all the strings operators 
applied to remove the excitations in the bulk of the chain, see 
Eq. (^. Top: The first pj is fused with the left MBS. The total 
weight is 4. It leads to a successful recovery of the encoded 
logical qubit. Bottom: The first 'ip is fused to the second ip. 
The total weight is 9. This leads to a faulty procedure and 
the logical qubit is lost. 


D. Coupling to thermal bath 

The total Hamiltonian for the wire and the thermal 
bath is the one considered in Ref. |32j 

H = HsYHb^Hsb- 


( 10 ) 
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Here we choose the system Hamiltonian Hs = ^ see 

Eq. (^. This choice ensures that all the errors originate 
purely from thermal fluctuations. Bravyi and Koenig 
have considered the opposite regime where errors are 
solely due to Hamiltonian imperfections with fij ^ 0 and 
t 7^ |A|.^ In their scenario, they showed that the life¬ 
time of the Majorana code increases logarithmically with 
L. As we will see, this is also true for our thermal-bath 
model at large temperatures. 

The bath Hamiltonian is bosonic and take the generic 
form 


hb = J2^j, ( 11 ) 

3 


where are local bosonic operators associated with 
fermionic site j. 


The last term Hsb in Eq. ( 10 ) stands for the bath-wire 
interaction that we assume to be parity conserving^ 


Hsb = - E ® “ 1 ) = “* E ® I2j-il2j ■ 

j j 

( 12 ) 

This form of the coupling seems quite natural since it 
corresponds to quantum fluctuations of the chemical po¬ 
tential. Note that 7/1 excitations are created in pairs by 
the bath since [727-1727, Sk] = 0 for /c 7^ 7,7 — 1 and 

{l 2 j-ll 2 j,Sj-i} = {^ 2 j -172 j,Sj} = 0 , ( 13 ) 

with 7 = and Sl = Sq = i^i^2L is the parity of 

the logical qubit. 


1. Davies Prescription 

Following the prescription of Davies that 

has become standard in many quantum information 
problemswe derive the master equation for the wire 
in the memoryless (Markov) limit, 

Ps{t) = -i[Hs. Ps{t)] + V{ps{t)). ( 14 ) 

The first term describes unitary evolution while the sec¬ 
ond one the exchange of energy between the bath and 
the wire. The so-called dissipator is 

V{ps{t)) = 7*(w)P5(i)(V(w))^ 

i,7 uj 

-i{(V(:^))tA7:^),ps(i)}) ,(15) 


The jump operators Ai{uj) cause transitions between 
eigenstates of i^s, with energy difference uj. We distin¬ 
guish between the following categories of transitions, see 

Fig.|^ 

• Pair creation (annihilation) of 1/1 in the bulk, with 
UJ = — 4 |A| {uj = + 4 |A|). 

• Pair creation (annihilation) of 1/1 at the boundary, 

with UJ = —2|A| (uj = +2|A|). More precisely, 
one 1/1 is created (annihilated) and the eigenvalue 
of 5 'o = changes sign. 

• Hopping of a 1/1 to a nearest-neighbor site inside the 
bulk, with cj = 0. 

• Hopping of a 1/1 into a nearest-neighbor MBS, with 
cj = +2|A|. 

• Hopping of a 1/1 out from an MBS to a nearest- 
neighbor site of the bulk, with uj = —2|A|. 

Here we use the convention that a negative sign of uj 
describes an energy transfer from the bath to the wire. 



FIG. 5. Pictorial box representation of all the possible error 
processes generated by the environment. Here we consider 
the simple time-independent problem where the MBSs are 
immobile at the end of the wire. On the right is shown the 
energy costs of the processes. A positive energy means that 
energy is transferred from the wire to the bath. 


where (uj) = (B\(s)Bj(tS)) are the bath 

spectral functions. Here, (• • •) = Tr('''e“^^^) is the 
thermal expectation value at inverse temperature [ 3 . 

In Appendix we derive explicit expressions for the 
jump operators A^(uj). Importantly, they are local and 
satisfy detailed balance. The Davies prescription en¬ 
sures that the steady state of Eq. ( 14 ) is the Gibbs state 

PGibbs = exp j /Tr (t 




The time evolution of the diagonal elements of Ps(t) 
decouple from the off-diagonal elements, see Appendix [Aj 
and the Pauli master equation for the population P(n^T) 
in an eigenstate |n) of Hs satisfies 

= W [W{n\m)P{m, r) - W{m\n)P{n, r)] , 
dr 

m 

( 16 ) 
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with transition rates 

W{n\m) = 'y{umn)\{'in\^'”"'i‘^mn)\n)f . (17) 

Here ujmn is the energy difference between the eigenstates 
\m) and \n). Note that \m) and \n) can be degenerate, 
with uJmn = 0- We have removed the superscripts on 
the spectral function as it does not depend on the 
position; we assume that the sites are coupled to identical 
and independent baths. 

In this work we consider an Ohmic bath where the 
rates j{uj) are 


After a time step r, we assume that a number n of 
spin flips has been applied on randomly chosen sites, 
where n is taken from a Poisson distribution with mean 
A^o = ^ Wtot • Here Wtot is the total rate of all al¬ 
lowed error processes. We assume that Wtot is state 
independent; this corresponds to an infinite tempera¬ 
ture scenario where 7 ( 0 ) = 7 (± 2 |A|) = 7 (± 4 |A|), see 
Eq. ( p!8| ) in the limit P ^ 0. For simplicity we choose 
7 ( 0 ) = 7 (± 2 |A|) = 7 (± 4 |A|) = Wtot/L =: ^tot- Since 
events in a Poisson process are independent^ one can sim¬ 
plify the problem by just considering a single spin. We 
have 




00 

1 — exp(—;dci;) 


(18) 


with coupling constant k and inverse temperature /3. 


(S^) = - = _g-2r«,tot _ (22) 

k=0 

Similarly, the standard deviation is given by 


E. Lifetime of Majorana Code: Infinite 
Temperature 

We focus here on the infinite-temperature limit, where 
we derive transparent analytical results for the lifetime 
of the Majorana code. As mentioned in Ref. [ 32 I the Ma¬ 
jorana code represents a useful quantum memory with 
a lifetime that grows with the wire’s size L. A similar 
scaling behavior was discovered by Bravyi and Koenig 
in Ref. [16] and Kay in Ref. |40l However, these refer¬ 
ences considered unitary evolution, while we focus here 
on dissipative dynamics. Unfortunately, the scaling is 
logarithmic and thus very modest. Here we present an 
analytical proof of this result. 


1. Analytic derivation 

It is convenient to map the problem to spins via a 
Jordan-Wigner transformation 

aj = ^ j S+ and 7 = 2a]aj - 1. (19) 

In spin language takes the simple form 

L-l 

= -4|A| ^ 7571. (20) 

i=i 

In spin language the logical states |0) and |1) are recog¬ 
nized as the states with all spins pointing along x and 
—X. The logical Z Pauli is then obtained by application 
of the parity operator 


L 

(21) 

i=i 

We model the error process taking place on the spin 
chain as follows. All spins point initially along —x. 


= •\/((7)2) - (7)2 = yi - e-4™tot . (23) 


We thus have 


l^tot = (7ot) = (E5f 
\i=l 

(24) 

From the central limit theorem we derive the probability 
distribution g{s) of the random variable s = 5'^°^ with 
standard deviation atot = VlVI — namely 



L) 


1 (s fJ'totA 


^tot ^TT 


(25) 


Error correction in the spin chain is performed 
straightforwardly: one flips either all the spins point¬ 
ing along —X or along x, such that all the spins point 
along the same direction after the error correcting step. 
The choice of which spins to flip is done according to 
the total number of spins that need to be addressed; we 
choose to flip the minimal number of spins. In the origi¬ 
nal fermionic language, this is the Minimal Weight Per¬ 
fect Matching algorithm of Sec. H C and the two choices 
correspond to the possibility to fuse the first -0 with the 
left MBS or not. As a direct consequence, the question 
whether error correction is successful after time r maps to 
the question whether the majority of spins still points in 
the same direction as the initial state. If, for example, the 
initial state has all spins pointing along the —x-direction, 
then the condition for successful recovery after time r is 
that (with high probability) 


L 

^7<0. (26) 

The times r for which error correction has a high proba¬ 
bility of success are thus those for which 


i^tot) + ^tot ^ 0 


(27) 
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and so 


L > - 1. (28) 

Stated in other terms, the Majorana code has no thresh¬ 
old; the lifetime of the memory increases only logarithmi¬ 
cally with L. A similar scaling behavior was discovered 
by Bravyi and Koenig in Ref. ITbl and Kay for the surface 
code with one-dimensional Hamiltonian perturbations.^^ 
Note that these references considered unitary evolution, 
while we focus here on dissipative dynamics. The infi¬ 
nite temperature limit is a special case which allows for 
straightforward analytic treatment. However, because it 
is in a sense a worst case scenario for error correction we 
expect the ” no threshold” result to hold for more general 
error models. 


Pfailure ) 



FIG. 6. Probability of failure Pfaiiure(T) as function of time for 
Kitaev wires of lengths L = 9,19, 59. The results are obtained 
in the high-temperature regime with Wtot = T|A|. The two 
MBSs are immobile and we perform the MWPM algorithm 
to remove ex cita tions. The solid lines represent the analytical 
result of Eq. (29) while the dots are obtained via Monte Carlo 
simulation. We find good agreement between analytics and 
numerics; the improvement becomes better for larger L, as 
expected. 


Finally, the probability that error correction fails after 
time r is then simply given by 

7'failure('r, T, L)ds = I + Erf ^ ^ 

''‘T29) 

In Sec. HE 2[ we compare our analytical result with a 
Monte Carlo simulation and find very good agreement. 


simulation by adding a -0 in boxes 2j — 1 and 2j. An 
even number of -0 in a box is identical to vacuum. Note 
that a 7/^ and an MBS can coincide in the same box. To 
be more precise, we implement the effect of the error 
operator 7172 ( 72 L-i 72 l) by adding a 7 /^ in boxes 1 and 
2 (L and T + 1), although box 1 (L + 1) carries an MBS. 
This is just to signify that the parity of the logical qubit 
has been flipped. So —Sq. A logical Z Pauli occurs 
when two -0 only are present in the chain, namely in boxes 
1 and 1/ + 1, see Fig. 

An iteration of the simulation decomposes into the fol¬ 
lowing steps, i) We register all the relevant parameters of 
the system, in particular the actual configuration of exci¬ 
tations. ii) For a given time interval 6t, if r + (5r ^ Tgim, 
we update the time to r Sr and go to step iii). If 
r Sr > Tsim we go directly to step v). The time Tgim is 
the simulation time and describes how long the wire and 
the thermal bath have been in contact with each other, 
iii) We draw the number n of error processes from a Pois¬ 
son distribution with mean Wtot Sr. It is worth pointing 
out that Wtot is a state-dependent quantity; for a given 
eigenstate |n) of Hs, the total transition rate is 

Wtot{n) = Y,W{m\n), (30) 

m 


where \m) are eigenstates of Hs. 

However, in the infinite-temperature limit considered 
here, the total error rate becomes state independent, iii) 
We apply n error processes randomly according to their 
relative rateP^ and go back to step i). v) We perform 
the MWPM algorithm described in Sec. m and finally 
record whether the error correction was successful or not. 

To obtain reliable statistics we perform these five steps 
on several thousands of samples for each Tsim- In Fig- 
we plot the probability of failure as function of time for 
different lengths of the wire. The solid lines describe 
the analytical results (29), while the dots are obtained 
from the Monte Carlo simulation just described. We see 
that both results coincide very well (with the agreement 
improving for bigger L) and the logarithmic lifetime of 
the memory is confirmed. 

At low temperatures (i.e. [3 1/|A|), the dominant 

processes leading to faulty error correction is the diffusion 
of a single pair of ^ particles. The reason is that, at low 
temperature, it is not favorable to create quasi-particle 
pairs and it costs much less energy for an existing pair to 
diffuse than for a new pair to be created. We will treat 
this case in the following sections, when we consider the 
trijunction. 


2. Monte Carlo Simulation for the wire 


III. TRIJUNCTION 


We have applied standard Monte Carlo methods to 
sample Eq. We use the box representation of Fig.|2] 
to describe the simulation. An error caused by a system 
operator 72 ^- 172^5 see Eq. (12), is implemented in the 


In this section we follow our earlier approaclP^ and 
present the main aspects of the trijunction setup. 

As a one-dimensional wire does not have enough space 
to exchange MBSs, Ref. [5] proposed to use a trijunction 















and to move MBSs by tuning locally the different chem¬ 
ical potentials. Reference [5] demonstrated that, when 
MBSs are exchanged in the trijunction, they obey the 
same non-abelian braiding statistics as Ising anyons. It is 
thus very important to understand their properties when 
coupled to a thermal environment. In particular, below 
we will determine how the induced thermal noise affects 
braiding in a nontrivial way. 

The trijunction Hamiltonian is taken to be time de¬ 
pendent, 

2L 2L-1 

i=i i=i 

2L-1 

+ y] (Aa^a^+i + A*a]^^aj) - t{a[^^aL+i + al,+iaL/2) 
i=i 

+{AaL/2aL+i + (31) 

The sites correspond to the horizontal wire and 

the sites L + 1,..., 21/ to the vertical wire, see Fig. 
The last two terms in i^trij(^) describe the coupling be¬ 
tween the vertical and horizontal wires of the trijunction. 
In terms of Major ana operators, the trijunction Hamil¬ 
tonian becomes 

. 2L 2L-1 

Htrijir) = + i|A| y] I2jl2j+1 

i=i i=i 

+i|A|7L72L+i, (32) 

where the last term describes the coupling at the tri¬ 
junction point. Note that the time-dependent chemical 
potentials satisfy /ij(r) ^ 0 at all times r. 

The set of chemical potentials fij (r) is controlled exter¬ 
nally such that the full braid depicted in Fig.[^ is imple¬ 
mented. In the following, topological segments are char¬ 
acterized by fij = 0 with Hamiltonian ^ see Eq. ([^. 
Nontopological segments have \fij \ |A|, such that they 

are well approximated by see Eq. 

Eor clarity, Eig.[^ shows the box representation of the 
full trijunction when the horizontal wire carries three 
MBSs i.e., in braiding stage Hi) of Eig. [^. Indeed, one 
must pay attention that a Majorana mode of the hori¬ 
zontal wire will be paired with a Majorana mode of the 
vertical wire during this braiding and this is reflected in 
the box representation as shown in Eig. 


A. Encoding 

Consider four MBSs 71 , 2 , 3,4 as in Eig. Eollowing 
the procedure of Ref. |49l we encode the logical qubit 
in a fixed-parity sector, say ^7172 ^7374 = +1- Thus, 
while the ground-state subspace is fourfold degenerate, 
we use only two states to encode the qubit. This is 
required since the overall parity is fixed by the super¬ 
conducting pairing terms, so that gate operations can 


“h 2 L 



FIG. 7. a) Pictorial representation of the trijunction setup. 
It is composed of a horizontal wire coupled to a vertical wire 
through hopping and superconducting pairing. The black 
(larger) dots represent the fermionic sites and the black bonds 
represent hopping and superconducting pairing. The red 
(smaller) dots describe the Majorana modes and the lines in- 
between their pairings. The four Majorana modes 71 , 2 , 3,4 are 
depicted, b) Representation of the braiding motion considered 
in this work. MBSs 2 and 4 remains immobile, while MBSs 
1 and 3 are exchanged. The black and gray regions represent 
respectively topological and non topological segments of the 
trijunction. 

only be performed within a fixed-parity sector. The log¬ 
ical qubit states satisfy 2717210 ) = 2737410 ) = | 0 ) and 
ni72|i) = *7374|i) = -|i)- 

Again, the logical X, T, and Z Pauli operators are rep¬ 
resented in terms of -0 particles inside MBSs, see Fig. 


IV. UNITARY EVOLUTION AND MBS 
MOTION IN THE TRIJUNCTION 

The motions of MBSs in a braiding sequence are per¬ 
formed unitarily. Therefore it is worth spending time 
to describe the unitary evolution of excited states when 
MBSs are moved. Indeed, it is essential to understand 
how moving MBSs interact with T/^-particles if one wants 
to simulate the dynamics of the system. The rules gov¬ 
erning the interactions between moving MBSs and excita¬ 
tions were reported in Ref. |32l here we present a detailed 
analysis leading to these rules. 

A. Adiabaticity 

For any braiding protocol to be valid, the MBSs must 
be moved sufficiently slowly with respect to the gap sep- 
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FIG. 9. Pictorial box representation of logical Z-, JX-, and Y- 
FIG. 8. a) Pictorial representation of the trijunction during errors in terms of '0-particles inside MBSs. 
braiding stage Hi) of Fig.[^; three MBSs lie on the horizontal 
wire, b) Box representation of a). 

one site to a nearest-neighbor site. 


arating the ground states and the rest of the spectrum; 
here this is the superconducting gap |A|. In other words, 
the chemical potential Hj{r) at site j must be varied 
slowly enough. This was a central assumption to our 
previous work Ref. isa here we give an explicit formula 
for the time-dependent chemical potentials. We imple¬ 
ment adiabaticity by choosing 

= -10"^|A|2t?(r|A| - Tj\A\) {t - Tj) , (33) 

where i^{r) is the Heaviside theta function and rj is fixed 
by the details of the braiding motions and determines 
when the chemical potential at site j starts to change. 
We have tested numerically whether the above functional 
form of the chemical potential is good enough to remain 
within the adiabatic regime. We have diagonalize a time- 
dependent four-site model; starting from the ground- 
state we have calculated its time evolution up to time 
r and its overlap with the instantaneous ground state at 
time r. The overlap was very close to 1 at any time; for 
example the value of the non-adiabatic matrix elements 
at time r|A| = 10^ was not larger than 1.6 x 10“^. 

We find that a good rule of thumb is that an MBS has 
moved from site j to site j + 1 when /ij = 10|A|; in other 
terms it takes a time r|A| « 10^ to move an MBS from 


B. Linear Motion 


For the linear motion along the horizontal or verti¬ 
cal wires, our analysis is based on numerical diagonal- 
ization of a small Kitaev wire composed of four sites 
where we successively decrease the chemical potentials 


on the different sites, see Fig. 10 r. This is done very 
slowly such that the system remains in an instantaneous 
energy eigenstate. We have computed numerically the 
time-ordered exponential 


U{0,r) = Texp ^—i J dr'Hsir') 


(34) 


that describes the unitary evolution under Hs{r). Here 
T is the time-ordering operator. 

Consider the initial configuration of '0 shown in 
Fig. [T^, i.e. one -0 inside the leftmost MBS and an¬ 
other one in the second box, such that d[di = 1. 
The initial state is thus an eigenvector with eigenval¬ 
ues dj^do = l^d\di = I,d 2 <i 2 = O^d^^ds = 0. At 
time r = 0, the parity of the logical qubit is given by 
2dldo — 1 = '^7 i78 = +1- After decreasing the chemi¬ 
cal potential according to /ii(r) = —10 “^|Apr until it 
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FIG. 10. Four-site Kitaev wire model, a) The black dots 
represent the sites of the Kitaev wire. The chemical poten¬ 
tials and /J^ 2 {r) on the first two sites are successively 

decreased according to = — 10 “^|A|^t and = 

—10“^|A|'d(T|A| — 10 ^)(t|A| — 10^). The solid black lines 
describe hopping and superconducting pairing. Below you 
can find the Majorana representation of the same wire, where 
smaller (red) dots represent Majorana modes. The double 
lines describe the varying chemical potentials, b) Box repre¬ 
sentation of the four-site model. The '0- and ^'-excitations 
are explicitly shown and gray shaded boxes correspond to 
nontopological segments. The initial state supports two -0- 
excitations that evolve during the unitary evolution i) and 
ii). 


reaches the value /ii(10^/|A|) = — 10|A|, the MBS has 
moved to the right and the parity of the logical qubit 
is changed to iysys- The amplitude of the chemical po¬ 
tential on site 1 being large, the operator 27172 becomes 
close to an eigenoperator of the Hamiltonian. 


In Fig. we plot the expectation values of 27172 and 
^7378 as function of time. As both go to +1 at time 
r|A| = 10 ^, we interpret the results as follows: the ip 
carried by the MBS stays bound to the MBS, while the 
other is transferred from the topological segment into 
the nontopological one. An excitation in a nontopologi¬ 
cal segment is called a pj' to notify that it has different 
attributes, e.g. a higher energy. In order to localize the 
'^'-excitations, the chemical potentials in the nontopolog¬ 
ical segments have a gradient as shown in Fig. [^. 


Let us now decrease the second chemical potential ac¬ 
cording to /22(^) = —10“^ ^ while, at the same time, 

we continue to decrease / 2 i(r) until time r|A| = 2 x 10^, 
see Fig. jl^. We thus have /ii(2x 10^/| A|) = —20 |A| and 
/i 2 ( 10 ^/[A|) = — 10|A|. The parity of the logical qubit 
becomes 27578. We expect to see the ip' excitation im¬ 
mobile in the nontopological segment, while the ip bound 
to the MBS moves together with the MBS further to the 
right. This is exactly what we observe in Fig. 12 where 


we have plotted (27172) and (27573) as function of time. 
We point out that it is necessary to maintain different 
chemical potentials on the different sites of the nontopo¬ 
logical segment in order to localize the '^'-particles, see 

Fig. [Tjj. 



FIG. 11. Expectations of 27172 and 27378 as a function of 
time for the four-site model of Fig. The plots correspond 
to transition 2 ) of Fig, p^. We see that both quantities goes 
to 1 as / 2 i(t) is decreased from 0 to — 10|A|. When the chem¬ 
ical potential on the first site is as low as — 10 |A|, then the 
MBS has moved to the right and the parity of the Majo¬ 
rana qubit becomes 27378 - The fact that ( 27172 ) ^ 1 and 
( 27378 ) ^ 1 show that the -0 bound to the MBS stays bound 
to the MBS while the other '0-excitation is transferred to the 
nontopological segment and thus becomes a '^'-particle. 

We have performed several similar tests, starting from 
different configurations of ' 0 -excitations. All the conclu¬ 
sions are the same and can be summarized in terms of the 
following rules, i) When an MBS moves into a topological 
segment and crosses a ' 0 -excitation, then the ' 0 -excitation 
is transferred to the first site to the left of the MBS into 
the nontopological segment and stays immobile, ii) A '0 
inside an MBS moves together with the MBS. In case 
of reverse motion, i.e. when the MBS moves into the 
nontopological segment, then a ' 0 ' from the nontopologi¬ 
cal segment will be transferred back into the topological 
segment. We have represented these two rules pictorially 
in Fig. It is important to recognize that while the 
overall parity of the system is conserved, the parity of in¬ 
dividual topological segments is not preserved. This is a 
crucial difference as compared to the previous case with 
two immobile MBSs, as now logical X and Y errors are 
possible. 

In Ref. [33] we summarize all the unitary evolutions nee- 
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V. ADIABATIC DAVIES EQUATION 


We must take the time-dependence of the trijunction 
Hamiltonian into account when writing down the master 
equation!^ In the adiabatic limit, it is correct to gener¬ 
alize Eq. (14) to 


hi uj{t) 

(35) 


where uj{r) are the time-dependent energy differences in 
the spectrum of i^trij('r). 

The populations follow an adiabatic Pauli master equa¬ 
tion 


dP{n{T),T) ^ ^ [W{n{T)\m{T))P{m{T),T) 

m{r) 

-W{m{T)\n{T))P{n{T),T)] , (36) 


FIG. 12. Expectations of 27172 and 27578 as a function of 
time for the four-site model of Fig. The plots correspond 
to transition 22) of Fig. lOb. Here the second chemical po¬ 
tential, /22 (t), starts to decrease while the first chemical po¬ 
tential, / 2 i(t), continues to decrease. When / 22 (t) is negative 
enough, then the MBS has moved to the nearest-neighbor site 
and the parity of the logical qubit becomes 27578. The fact 
that (27172) = ! confirms that i/^'-particle remains bound at its 
position in the nontopological segment. Since (27573) ^ 1 as 
the chemical potential is decreased, we conclude again that 
the i/;-particle inside the MBS remains inside the MBS during 
the unitary evolution. 


i) 


• 



• 



• 


• 



FIG. 13 . Box representation of the unitary evolution of 
moving MBSs in the presence of excitations, i) When a 
MBS moves over an existing i/;-excitation, the i/^-excitation 
is transferred to the nontopological segment and becomes a 
i/;bexcitation. ii) A ' 0 -excitation trapped inside a MBS re¬ 
mains trapped during the motion of the MBS. 


essary to simulate the system; the unitary motions of an 
MBS over the trijunction point need to be obtained by 
simulating a minimal six-site trijunction model. 


with 


W{n{T)\m{T)) =7(a;m„(r)) |(m(T)|A*™"(a;m„(r))|n(T))| 

(37) 

The bulk error processes and the associated rates re¬ 
main the same as in the time-independent scenario, see 


Sec. HD 1 To be more precise, the error processes away 
from the moving MBSs, including at other MBSs that 
are for the time being stationary, are the ones presented 
in Sec. EDT] However, more complicated boundary pro¬ 
cesses appear because of the motion of the MBSs. In 
Appendix we present some examples. An exhaustive 
table of the more than 200 distinct allowed processes can 
be found in Ref. |33l 

It is worth pointing out that the system-bath interac¬ 
tion of Eq. (12) does not support the creation of excita¬ 


tions in the nontopological segments of the trijunction. 
Indeed, when the chemical potential is very negative, 
the eigenstates of a nontopological segment approaches 
the eigenstates of and the coincidence becomes 

better as the chemical potential becomes more negative. 
Since ^ Hs b] = 0, creation of excitations in the 

nontopological segment is suppressed as the chemical po¬ 
tential decreases. However, this does not mean that no 
excitations will ever be present in the nontopological seg¬ 
ments, as we discussed in Sec. |IV[ 

The time dependence of Hs{r) must also be taken into 
account to calculate the rates of all the error processes. 
Eor example, when changing the chemical potential from 
time 0 to time r with r|A| ^ 10^, such that an MBS has 
moved by one site, one obtains the rates associated with 
the possible error processes by integrating Eq. (37); we 
defer a detailed discussion to Appendix |B| 
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VI. ERROR CORRECTION AND DANGEROUS 
PROCESSES IN THE TRIJUNCTION 

The error correcting procedure applied here is an adap¬ 
tation of the greedy algorithm proposed by Wootton in 
Ref. |50l We point out that this is not a MWPM algo¬ 
rithm in the sense described in Sec. cmi We have cho¬ 
sen this decoder because of its great simplicity. Also, 
choosing any other decoding scheme would not change 
the main message of our paper, as we will see. More 
generally, many matching procedures might be applied 
to the same model, each of which potentially having a 
different threshold.^ 

Our error-correcting scheme here is passive, meaning 
that we apply it at the end of the quantum computing 
protocol (that below will consist in a single braid only). 
This is in contrast with active error correction where er¬ 
ror correcting steps are performed in the midst of a braid 
sequence of a full quantum computing protocol. We be¬ 
lieve that it is useful to make the distinction between the 
following active scenarios: 

1. Error correction is executed at the end of each 
braid. 

2. Error correction is performed repeatedly during 
each braid. 

As we will see, our results imply that passive error cor¬ 
rection and even active error correction that is performed 
at the completion of braids, do not lead to a lifetime that 
increases with the size of the trijunction; a scheme where 
error correction is executed during each braid is required 
to cure the dangerous errors. 

Below we summarize the main steps of the error cor¬ 
rection algorithm for the trijunction: 

1. Loop through all sites of the trijunction to find the 
pairs of quasi-particles (t/^, or MBS ) that are 
at minimal distance k. Start with k = 1. 

2. If there a no multiple possibilities, annihilate the 
corresponding pairs. If there are multiple possibil¬ 
ities, e.g. in the case that three quasi-particles are 
positioned such that one of them is at distance k 
from the two others, we apply the following rules: 

• If all the excitations are on the horizontal wire, 
we pair the particles from left to right. 

• If all the excitations are on the vertical wire, 
we pair the particles from top to bottom. 

• If two excitations are on the horizontal wire 
and a third is on the vertical wire, we annihi¬ 
late the pair composed of the leftmost quasi¬ 
particle on the horizontal wire and the upper¬ 
most quasi-particle on the vertical wire. 

• If two excitations are on the vertical wire and 
a third is on the horizontal wire, we annihilate 


the pair composed of the uppermost quasi¬ 
particles on the vertical wire and the quasi¬ 
particle on the horizontal wire. 

3. If there are still some excitations in the bulk of the 
trijunction, repeat the procedure with /c + 1. 

To illustrate this procedure we present a pictorial repre¬ 
sentation of one error correction step in Eig. 



FIG. 14 . First iteration {k = 1 ) of the greedy decoder used in 
our trijunction simulations. After having identified the pairs 
of excitations at distance /c, we annihilate them according to 
the protocol described in Sec. im Our algorithm is a straight¬ 
forward adaptation of the one presented in Ref. 1501 


A. Dangerous Errors 

This section contains the central result of our work: the 
identification of so-called dangerous errors that prevent 
the lifetime of the Majorana trijunction from increasing 
with the system size. 

Consider the situation in which an MBS is moving and 
a pair of excitation is created, one inside the MBS and 
one inside the bulk of the trijunction, see Fig. While 
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the two are originally created as a pair in neighbor¬ 
ing boxes, the motion of the MBS drags along one of 
the ^ and separates it from its partner. In other terms, 
the braiding renders an originally local error source com¬ 
pletely nonlocal. We thus call an error process that cre¬ 
ates a 7p inside a mobile MBS dangerous. The effective 
non-locality of the noise prevents our algorithm from suc¬ 
cessfully recovering the stored quantum information. For 
example in Fig. [15^, a single error event will not be cured 
by our algorithm and will lead to a X error. Note that in 
Fig. we have drawn two error processes: a dangerous 
error process, where a is created inside an MBS, and 
an inoffensive error process where a pair of excitations is 
generated in the bulk of the vertical wire. 

A natural question that arises is whether a better al¬ 
gorithm could take into account the nonlocality of the 
noise in a clever manner. Unfortunately this is impossi¬ 
ble if error correction is not performed during braiding. 
The reason is that different error processes can lead to 
exactly the same error syndrome. In Fig. \TE^ and b, we 
depict two error processes that generate the same syn¬ 
drome. The main difference between them is the occur¬ 
rence of V^-particles inside MBSs. This can be traced 
back to the moment where a dangerous error happens. 
In Fig. it happens at the beginning of the braid, 
while it happens at the end of the braid in Fig. [T^. If 
one syndrome is successfully cured by an algorithm, the 
other one will lead to failure. A central ingredient for 
the emergence of such ambiguity is that different MBSs 
travel over the same segments of the trijunction during 
braiding, making it impossible to identify which bulk 
should be paired with which MBS. Since the probability 
of dangerous events is finite and does not depend on the 
size of the trijunction, we expect the lifetime of coherence 
of the trijunction qubit to be independent of L. 

We point out that the physics of dangerous error pro¬ 
cesses is the same at high 11 /^ 1 ) low ^ 

I/|A|) temperatures. Therefore, we expect the restric¬ 
tion due to dangerous error processes to be qualitatively 
identical in both regimes. 


VII. MONTE CARLO SIMULATION OF THE 
TRIJUNCTION 


To confirm our predictions, we perform a standard 
Monte Carlo simulation for the trijunction and determine 
the evolution of the stored quantum information under 
the adiabatic master equation ( [^ . We focus here on the 
low-temperature regime with P = 4/|A|. 

The Monte Carlo simulation consists essentially of 
the same five steps as in the simulation described in 
Sec. |IIE 2 | for the single Majorana wire. However, some 
care has to be taken in the low-temperature regime. In¬ 
deed, in that case the total probability Wtot that an er¬ 
ror event occurs, see Eq. (30), is strongly state-dependent 


and we cannot always approximate it by a constant. This 
is the case because the spectral function 7 ( 0 ;) depends 



FIG. 15 . Pictorial representation of two error sequences lead¬ 
ing to the same error syndrome. The difference between a) 
and h) is the occurrence of p inside MBSs. Since p trapped 
inside an MBS is invisible to an external observer, both sit¬ 
uations have the same error syndrome. Independently of the 
details of the error-correcting algorithm, if one situation is 
successfully corrected, the other one will lead to failure. 


very much on the value of uj at low temperatures. 

One possible way to solve this issue would be to apply 
an alternative set of five Monte Carlo steps: i) Register 
all the relevant parameters of the system, in particular 
the actual configuration of excitations, ii) Calculate the 
time 5t for the next error process to occur, drawing 5t 
from an exponential distribution oc exp(I/IFtot)- hi) Up¬ 
date the time to r + 5 t. li t ^ 5 t ^ Tsim, go to step 
iv). Otherwise go directly to step v). iv) Apply an er¬ 
ror event randomly according to their relative rates. Go 
back to step i). v) Perform the error corecting algorithm 
described in Sec. |VI| and finally record whether the error 
correction was successful or not. 

Such a procedure is perfectly valid in the low- 
temperature regime, but only when the MBSs are im¬ 
mobile. Indeed, when MBSs are in the process of being 
moved this method cannot be applied. The main problem 
resides in the fluctuating 5t drawn from the exponential 
distribution. When MBSs are moved, one needs to de¬ 
fine a time-step Stm for an MBS to be carried to the 
nearest-neighbor site. For example, here we have chosen 
(5rM|A| = 10^. It is however clear that, most of the time, 
Sr drawn from the exponential distribution would never 
be an exact multiple of Stm and thus we cannot decide 
by how many steps the MBS must be moved during the 
time interval Sr. Therefore, this method is applicable 
only when the MBS are kept immobile; so, this method 
is applicable during the time in between braids. 
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During the braidin g moti ons, we apply the sequence of 
steps i)-v) from Sec. HE2 but by taking care that VEtot 
is state dependent and by ensuring that the number of 
errors n drawn from the Poisson distribution is either 0 
or 1; we assure this by choosing a small enough coupling 
constant n = 2 x 10“^|A|. We point out that this choice 
of n is only relevant for our numerical procedure to be 
valid, but it does not hide any important physical issue. 
In fact, allowing n > 1 would be problematic (in the low- 
temperature regime) since each time an error is applied, 
the change in total probability Wtot is very drastic and 
must be taken into account. 

In Fig.[^ we present the probabilities px^ Py^ and pz 
that a logical X—, F—, or Z—error occurs for trijunctions 
of various lengths, see also Ref. |32j When we assume 
perfect error correction, the noise acting on the logical 
qubit is unital, i.e. the fixed point is a completely mixed 
state, and takes the form of a generalized depolarizing 
channel^ 


Ptrij ^ ^ ^ PO O Ptrij O , (^8) 

0=PX,Y,Z 

where / is the identity operator and pi -\-px -\-PY -\-pz = 1 • 
This is in contrast with the underlying physical noise that 
is not unital here. 

We observe that px and py significantly increase dur¬ 
ing the braiding time, while pz remains small. This is 
due to the presence of dangerous error processes during 
braiding as explained in Sec. |VI A[ It is also interesting 
to distinguish between X and Y errors. It is clear from 
the plots that px increases faster than py and there is a 
period of time where py does not increase. In fact, con¬ 
sidering only single-error events (a very accurate approx¬ 
imation at short times) the environment cannot produce 
a Y error during the first L/2-\-l time steps. Indeed, the 
only possible configuration of T/^-excitations correspond¬ 
ing to a F-error after error correction, resulting from 
a single error event during braiding, is the one in the 
bottom left of Fig. Such a situation arises when a 
dangerous error happens during the braiding stage i) of 
Fig-lZl However, the V^-particle created in the bulk of 
the horizontal wire must he closer to the right boundary 
than to the left boundary. Therefore, such a F-error in¬ 
deed cannot happen during the first I//2 +1 time steps of 
braiding. On the contrary, an X error can be produced 
by a single error events at any time during the braid. We 
also point out that dangerous error processes are more 
probable than creation of a pair somewhere in the bulk, 
since the energy cost is lower. A transition energy of 2|A| 
vs. 4|A| makes a considerable difference, and the strong 
dependence of on uj also plays a role. 

In Fig. [^we plot px-\-Y = Px Fpv for different L. Af¬ 
ter the end of the braid px-\-Y remains constant since no 
dangerous errors are possible anymore when MBSs stay 
immobile, although X and F errors can be interconverted 
during this period. Indeed, X and F errors are possible 
only when the parity of topological segments is broken 
and this is solely possible when MBSs move. The most 


important feature of the growth of Px+Y is that it is 
completely independent of L at small times; the lifetime 
of the memory does not grow with L. This is in com¬ 
plete agreement with our discussion of dangerous errors 
in Sec. \VTA\ It is worth pointing out that the braid¬ 
ing time grows linearly with the size of the trijunction. 
Therefore, the probability that a dangerous error occurs 
during braiding is higher for a larger trijunction. This is 
observed in Fig. 17, where Px+Y at the end of the braid 
is bigger for larger trijunctions. 

At small temperature, the origin of a non vanishing pz 
probability is the creation of a pair of t/^’s that diffuse 
across the trijunction. However, it takes more time for 
a pair to reach the MBSs when the trijunction is longer, 
therefore we expect pz to decrease with increasing L, and 
this is the case in Fig. 17 r. 

In the context of a full quantum computing protocol, 
where several braids are executed, our results show that 
performing error correction either at the end of all the 
braids or at the end of each braid is not enough to cure 
the failure induced by dangerous errors. 


VIII. CONCLUSIONS 

In this work we have investigated the self-correcting 
properties of Majorana ID quantum computing architec¬ 
tures. In particular, starting from a microscopic model, 
we focused on the situation where MBSs are braided in a 
trijunction setup coupled to a parity-preserving bosonic 
environment. 

While a single wire with immobile MBSs represents 
a truly self-correcting quantum memory with a lifetime 
that increases with the size of the wire, this is not true 
anymore when MBSs are exchanged in a trijunction ar¬ 
chitecture. The main reason is the occurrence of so-called 
dangerous errors that are solely due to the motion of 
MBSs; an MBS can trap an excitation and drag it along 
during braiding, thus rendering a local source of noise 
highly nonlocal in its effect. In this case, error correction 
at the end of the braid is insufficient to recover the stored 
quantum information. 

In the context of a full quantum computing protocol, 
where several braids are executed, our results imply that 
passive error correction (at the end of all braids) and even 
active error correction, in which correction is performed 
at the end of each braid, is too weak to counteract the 
negative effects of dangerous errors. The only possibil¬ 
ity we envision to preserve the stored quantum informa¬ 
tion is to perform repeatedly error correction during each 
braid: in the very simple example of Fig. performing 
error correction both at stages Hi) and v) would allow 
successful recovery of the logical qubit. Our results are 
in agreement with the more and more popular view that 
active error correction is necessary even in non-abelian 
topological systems .1^®^ In light of this discussion, we 
point out that increasing L, far from improving the life¬ 
time of the trijunction logical qubit, actually makes the 
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FIG. 16. Plots of px, Py^ and pz as function of time. MBSs 
1 and 3 are exchanged while MBSs 2 and 4 remain immobile, 
see Fig. The time to braid is proportional to the size 
of the trijunction and we have Tbraid = (2L + 2)10^/|A|. For 
T < Tbraid, i-G. wheu the exchange is not finished, we exchange 
the MBSs until time r with the coupling to the thermal bath 
being on. Then we unitarily (no coupling to the thermal bath) 
finish the braid and perform error correction at the end. Here 
we use — 2 X 10“^|A| and 13 = 4/|A|. 



FIG. 17. a) Probability px+Y (solid) and pz as function of 
time for trijunctions of length L = 19,29,49. The details of 
the plot are the same as in Fig. |16[ b) Artistic representation 
of the probability of failure as function of time when three 
consecutive braids are executed. The time interval between 
two braids is Sr and the braiding time is Tbraid- 

the same restrictions as the ID case studied here. In 
2D setups, the paths followed by the braided MBSs must 
cross at least once, but do not need to overlap over a large 
region. Therefore, dangerous error processes should oc¬ 
cur less frequently than in ID implementations; but we 
keep in mind that a single uncorrectable error is enough 
to prevent the lifetime of the topological qubit from in¬ 
creasing with the system size. 
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situation worse because the time to braid MBSs adiabat- 
ically is proportional to 

We also comment that while our results put restrictions 
on the self-correcting properties of this specific quantum 
computing architecture, one can expect ot her sc hemes, 
such as interaction-based braiding of MBSs,l^^l^ to be¬ 
have in a more favorable way. Einally, it is a priori not 
clear whether braiding MBSs in 2D systems suffers from 
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Appendix A: Davies Prescription: 
Time-Independent Hamiltonian 


In this Appendix we aim to derive explicit expres¬ 
sions for the jump operators appearing in the mas¬ 

ter equation (15) for the time-independent problem; for 
simplicity we focus here on a single wire. The time- 
dependent case is treated in an exactly similar fashion 
because we consider only adiabatic time variation. 

Let us start from the system bath Hamiltonian 


Hsb = -2^Bj 




(Al) 


where the constant in Eq. (12) has been ignored because 
it only leads to a renormalization of the bath Hamiltonian 
Hb‘ Rewrite Hsb in terms of the eigenoperators dj and 
dj that diagonalize see Eq. We have 

Hsb = -i2Bi (g) (do + d\){di + d\) 

\L/2\-l 

— 2 B2j (g) {d2j-l — + <^2j) 

i=i 

rL/21-1 

-2 X ® id2j - dlj){d2j+i + 4j+i) 

i=i 

0 (di-i - 4-i)((^o - 4) • (A2) 

It is now useful to decompose Hsb into three physically 
relevant terms: 

Hopping: 

rL/21-i 

^hopping •— i2B\ (g) d()d^ 2 ^ ^ -^2jf C) d2j — '\_d 2 j 

i=i 

rL/ 21-1 

-2 ^ B 2 j +1 0 C^2j4i+1 “ 0 dL-idl + h.c. 

j=i 

(A3) 


Pair creation: 


A 


creation • — 


rL/ 21-1 

-i2Bi 0 dld[ +2 ^2 0 4i-i4i + B2j+i 0 <^2/4/+!) 

i=i 


+i2BL 0 4 - i 4 

(A4) 


Pair annihilation 


^annihilation • — 


rL/21-1 

i2Bi (g) dido H- 2 E {B2j (g) d2jd2j-i H- ^2^+1 C) d2j+id2j) 

i=i 


-i2BL (g) dodp-i 


(A5) 
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Following the Davies prescription, we calculate the 
Fourier transforms of the above operators by first writing 
their time evolution with respect to Hs^ 

^iHst g-iffst 

= \m){m\e'‘^^^\k){k\A^\£){£\e~‘‘^^^\n){n\ 

=, (A6) 

m,n 


Creation terms: 


^Lation(-2|A|) = 

E 

\m){m\2d\d\\n){n \, 


Cm—Cn = — 2 A 


AcGation(-2|A|) = 

E 

\rn){m\2dld\_^\n){n\, 


Cm—Cn = — 2 A 


4tation(-4|A|) = 

E 

\m){rn\2d\^_^d\^\n){n 


e^-e^ = -4|A| 


4ttkn(-4|A|) = 

E 

\rn){rn\2d\p\^j^-^\n){n 


Cm 



Annihilation terms: 


where \i) is an eigenbasis of Hs with eigenenergies and 
C G {hopping, creation, annihilation}. We have used the 
decomposition of unity 1 = "^^\Tn){m\. The Fourier 
components of A(^ are then simply given by 


^c(^)= \m){m\AQ\n){n\. (A7) 

Cm Cn 


^annihilation I ^ I) E |m)(m|2di(io|n)(n| , 

Cm Cr!,=2|Z\| 

^annihlation(^l^l) E \m){m\2dL-ido\n){n\, 

Cm Cn=2|Z\| 

^annihilation(^l^l) = ^ \m){m\ 2 d 2 jd 2 j-i\n){n\, 

€„^-e„=4\A\ 

^anntilation(4|A|) = \m) {m\2d2j+ld2j\n) {n\, 

e^-e^=4|A| 

(AlO) 


We can then easily identify the following relevant system 
operators: 

Hopping terms: 


■hopping(“2 A ) — 

m)(m 2(iod| n)(n , 


Cm—Cn = —2| A| 

■hopping(-2|A|) = 

y] \m){m\2dod\_^\n){n 


Cm ^n— 2|A| 

-hopping(2|A|) = 

\m){m\2didl\n){n\ , 


Cm Cn—2|A| 


with j = 1,..., \L/2] - 1. 

In the box representation, see Fig.[^ the terms with en¬ 
ergy argument ±2|A| correspond to hopping at the ends 
of the wire, where the 'ip hops out from the MBS to inside 
the neighboring box or vice versa, or to a process at the 
boundaries where a is created or annihilated inside the 
MBS in conjunction with a second at the neighboring 
box inside the bulk of the wire. The terms with energy 
argument 0 are associated with hopping processes inside 
the bulk of the wire, where a hops from one box into 
another one without energy cost. Finally, the terms with 
energy argument ±4|A| correspond to processes where a 
pair of excitations is created or annihilated in the bulk 
of the wire. 


1. Pauli Master Equation 


^hopping(2 |A|) = \rn){rn\2dL-idl\n){n\, 

Cm Cn=2|A| 

^hopping(O) = E \m){m\2d2j-id\^ + 2d2jd\^_^\n){n \, 

e^-e^=0 


4oppkg(o) = E i”i)(”i|2(^2i4i+i + 2d2j+id2_^-|n)(n| , 

e^-e^=0 

(A8) 


We present a proof that the diagonal terms decouple 
from the off-diagonal terms in the master equation ( p^ , 
leading to the Pauli master equation (16) for populations. 

Consider the eigenstates |n^) of such that Hs\ni) = 
en\ni). Here i = 1, 2 ,... indexes the degeneracy of the 
energy level. 

Assumption: If {ma\A\^{uj)\nk) 7 ^ 0, (here Cm — = 


uj) then there are no other system operators A\^, that 
cause a transition between \nk) and any of the degenerate 
states \mj 3 ) with energy e^. 

Note that it is straightforward to see that the above 
assumption is satisfied in our case. We thus have. 


with j = 1 ,..., \L/2'\ - 1. 
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inul'DiPsMnk) = Y1 ({nk\Ai{oj)psAl,,{co)^nk) - l{n,\Al^,{co)^Ai{aj)ps\n,) -\{nu\psAi^,{u:)^A^^{u:)\nu)\ 

i,j,r],r]' uj ^ ' 

= ^ 7*^H I E (?lfc|^Uw)|m„)(ma|/)s|K;3)(K/3|^^/MVfc) - E I 

uj,i,j,r],r]' \rn,K,a,l3 j 

E EE^'^'(^)(«'^i'^^i"^)("^i4HV«)(rd4(^)i«o (All) 

77,77' r,s (5,^ 

= E 7“(w)(’^fc|A(,(w)|m„(i,r?))(TOa(i,r?)|p5|TOa(i,r?))(m„(i,r?)|^(,(w)dnfc) 

01,7,77 

E 7”(w)(nfc|(A^(w)dma(z,r/))(m„(i,7?)|A^(w)|nfe)(nfe|p5|nfc) 

01,7,77 

-i ^ 7”(w)(nfc|/9s|nfc)(nfe|^),(w)dma(i,7?))(m„(z,r/)|^(,(a;)|nfc). (A12) 

01,7,77 


The Pauli master equation (16) follows then directly. 

Here we have used the orthonormality relation 
(^al^/ 3 ) = ^mn^ai3- The State \mo,) in the sums ( [All ) 
that has non vanishing matrix element {mcx\A\^{uj)\nk) 


depends on the system operator A\^{uj). For the sake 
of clarity we have thus introduced the notation IrUc^) = 
\ma{i,p)). 


Appendix B: Error Processes: Time-dependent Case 

The calculation of the rates of the error processes hap¬ 
pening during the motion of MBSs must be performed 
with great care. In the following we present a detailed 
analysis of two cases. All the remaining cases shown in 
Ref. [ 33 ] are treated similarly. 


1. Linear case 


The linear case is easily understood in terms of a four- 
site model where the chemical potential on the first site 
is varied according to /ii(r) = — 10“^|Apr. The red dots 
in Fig. [ 1 ^ represent the Majorana modes, while the sin¬ 
gle solid lines represent the coupling between them in 
the Hamiltonian. The double solid lines represent the 
time-dependent chemical potential. We assume that a 
region is nontopological when the chemical potential sat¬ 
isfies |/i| > 10|A|. Therefore, we say that the MBS has 
moved to a nearest-neighbor site after a time r| A| > 10^. 
The black arrows between different branches of the spec¬ 
trum in Fig. \TEjp describe transitions caused by the sys¬ 
tem operators. These transitions are also shown in box 
representation in Fig. [T^. 

The time dependence of the problem requires that the 
rates are obtained by integrating Eq. (37) over time. The 


upper integration bound T must be big enough such that 
the MBS has moved, i.e. as mentioned above we choose 


a) 0 =® —& ©—© @—e © 

/Ui(r) = -10“^|Apr 



FIG. 18 . a) Pictorial representation of the four-site model; 
the dots represent the Majorana modes and the single lines in 
between represent the pairings. The double lines depict the 
time-dependent chemical potential on the first site. It varies 
according to — — 10“^|A|^t. h) Energy levels Ek{r) as 

function of time. The arrows describe the transitions above 
the vacuum produced by the system operator 7172, inducing 
transition z), and 7374 inducing transitions ii) and Hi), c) 
Box representation of transitions z), zz), and zzz). When the 
chemical potential is sufficiently negative, then the MBS has 
moved to its nearest-neighbor site and a pair of excitations 
has been produced. 


T|A| = 10^. We have 

W(n|TO)'"* = l J (iT|(m(r)|4i,2|n(r)|)7(wTO„(T)). 

(Bl) 
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We calculate numerically the integral (B1) by discretiz¬ 


ing the interval [0,T] and by replacing the integral by a 


to the linear case, the rates are obtained by integrating 
Eq. 


2. Trijunction 


The processes that happen at the trijunction are ana- 
lyzable with a six-site model, see Fig. 19 r). In Fig. 19 3 we 
show the spectrum as function of time. The black arrows 
identify three transitions caused by the system opera¬ 
tor 7374 , see the box representation in Fig. [T^. Similar 



FIG. 19. a) Six-site trijunction. The dots represent Ma- 
jorana modes and the single lines the interactions between 
them. The shaded area shows that a Majorana mode is de¬ 
localized over the trjinaction. The double lines represent 
the time-dependent chemical potential varying according to 
//(r) = — 10 “^|A|^t. b) Spectrum as function of time. The 
black arrows represent the transitions above the vacuum in¬ 
duced by the system operator 7374 . c) Box representation 
of the two transitions i) and ii). The red shape on the left 
trijunction represents the delocalized MBS at the trijunction. 
Note that the two Majorana modes on the vertical wire are 
localized on neighboring boxes because we diagonalize here a 
small trijunction. In a longer trijunction, however, the third 
box on the vertical wire would correspond to a fermionic mode 
and the fourth MBSs would lie at the bottom of the vertical 
wire. In order to calculate all the rates, it is enough to diag¬ 
onalize such a small trijunction but one must keep in mind 
that creating a '0 inside the MBS at the bottom box of the 
vertical wire would correspond to creating a bulk ' 0 -excitation 
in the longer trijunction (and thus would imply some cost of 
energy). 





























Supplementary Material to ^^Monte Carlo studies of the properties of the Majorana 
quantum error correction code: is self-correction possible during braiding?” 
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The goal of this supplementary material is to present an exhaustive list of the hundreds of processes that are 
relevant for the Monte Carlo simulation of the trijunction Hamiltonian in Eq. (30) of the main text. We use our box 
representation to describe all the relevant processes. We present both the unitary evolutions as well as the dissipative 
error processes that are caused by the system-bath coupling Hsb^ see Eq. ( 12 ) of the main text. Of course, we present 
only the nontrivial unitary evolutions while more unitary evolutions, where a Majorana Bound State (MBS) is moving 
and excitations are inert, have to be taken into account in the simulations. Compared with the case where MBSs are 
immobile, many more processes must be taken into account when the MBSs are braided. 

Eor the sake of clarity, we explain the language we will use below: when we say that an MBS is moving over the 
trijunction point, we have in mind the moment of braiding corresponding to Eig. 1. 





FIG. 1. Pictorial representation of the trijunction setup supporting four MBSs 71,..., 4. When we say that an MBS is moving 
over the trijunction, we have in mind the precise moment of the braiding where 71 moves from the horizontal wire into the 
vertical wire. 

We distinguish among the following cases where an error process can happen: 

1 . Inside a topological segment, away from MBSs. 

2 . Close to immobile MBSs. 

3. Close to mobile MBSs that are moving in the direction of the topological segment. In other words, the considered 
topological segment becomes shorter. 

4. Close to mobile MBSs that are moving in the direction of the nontopological segment. In other words, the 
considered nontopological segment becomes shorter. 

5. Close to a mobile MBS during its motion over the trijunction point. 

6 . Close to a mobile MBS just before its motion over the trijunction point 

In the following tables, the cases enumerated above i.e. 1,2,3,4,5, and 6 , are written in the upper left corner. On 
the top of the table we also mention whether the depicted processes correspond to unitary evolution or to dissipative 
dynamics. Eor each case, the first figure, i.e. the one we call small models is a model of the system at the considered 
time during braiding. We use this model to identify and describe the system operators 7 ^ 7 ^ that are at the origin of 
the diffusive processes; we write on the left of each diffusive process the system operator causing the transition. 

Also, we would like to comment on the meaning of Eigs. 2 a and b, which illustrate elements occurring in the 
tables. The red shapes at the trijunction point represent a delocalized Majorana mode. Since the trijunction region 
is composed of two sites (site L/2 on the horizontal wire and site L + 1 on the vertical wire), it can support two 




2 


fermionic modes. Therefore, it is possible that the trijunction hosts a delocalized Majorana and a T/^-excitation at the 
same time, as depicted in Fig. 2b. It is worth pointing out that the Majorana modes contributing to a delocalized 
MBS cannot be nearest neighbors. For example, the delocalized MBS at the trijunction on page 15 has contributions 
from Majorana modes 5 and 9 but not from mode 4. 

We note that some of the figures in the tables below contain a single excitation or p;'). However, since the bath 
is parity-preserving, the total number of excitations cannot be odd. Therefore, when a single excitation is present, 
one needs to keep in mind that at least another inert undrawn excitation must also be present somewhere in the 
trijunction. 

Below the pj— and t/;'— particles are represented in a slightly different manner as compared to the main text. 
However, this difference is purely graphical and there is no difference at all between the quasi-particles represented 
below and the ones from the main text. 

a) b) 




FIG. 2. Pictorial representation of the moment in braiding when the MBS is moving over the trijunction point. The red 
shapes at the trijunction represent the delocalized MBS. While in a) the trijunction point carries only an MBS mode, in h) an 
additional '0-excitation is present. 
































































































































































































































































Case 3. 

Four-site model 
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Case 3. 

Unitary Evolution 
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Case 3. 

Error processes supported by the bath 
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Case 4. 


Unitary evolution 
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15 


We point out that the following trijunction processes only occur during stage iii) of Fig. 7b in the main text (or 
Fig. 1 of this supplement). In fact, when the MBS is moving over the trijunction during stage ii) of Fig. 7b in the 
main text, the error processes are exactly the same as the one occurring in a single (bent) wire. 

We recall that the delocalized MBS at the trijunction has contributions only from Majorana modes 5 and 9. 
Therefore, the system operator 7374 that we will consider below cannot create or annihilate a ^ inside the delocalized 
Majorana at the trijunction point; thus we only depict the case without a ^ inside the Majorana at the trijunction 
when we study the effect of 7374 . 
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Case 5. 


Error processes supported by the bath 
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